Small-world behavior in time-varying graphs 
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Connections in complex networks are inherently fluctuating over time and exhibit more dimen- 
sionality than analysis based on standard static graph measures can capture. Here, we introduce 
the concepts of temporal paths and distance in time-varying graphs. We define as temporal small 
world a time-varying graph in which the links are highly clustered in time, yet the nodes are at 
small average temporal distances. We explore the small-world behavior in synthetic time-varying 
networks of mobile agents, and in real social and biological time- varying systems. 

PACS numbers: 89.75.-k, 89.75.Hc, 89.75.Fb, 87.19.1j 
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In the last decade, the study of complex networks has 
attracted a lot of attention in the scientific community 
as various social, biological and technological systems 
can be represented and analysed as graphs [l|. Typi- 
cally, such systems are inherently dynamic, with the links 
changing and fluctuating over time. Human contacts or 
relationships change over time because individuals lose 
old acquaintances, acquire new ones, or move over geo- 
graphic space 2, 3]. Communication in man-made net- 
works, such as machine connections and social interac- 
tions over the Internet, takes place at specific points in 
time 0-1(1 . New links appear while some others disappear 
in the World- Wide- Web 7], in patterns of interactions 
among gene from microarray time experiments 0, Q or 
in functional brain networks [T3, [H| • The time evolution 
of a network by the addition, as well as the deletion, of 
nodes and links has been extensively modelled with the 
main purpose of reproducing, asymptotically, statistical 
properties such as scale- free degree distributions P, [l2| . 
With only a few notable exceptions (l3H15) . less focus 
has been given to the characterization of the dynamics 
of complex networks in stationary conditions. 

In particular, the small- world phenomenon, i.e., the 
fact that real networks have high clustering coefficient, 
while the typical distance between their nodes is small as 
in random graphs, has been investigated in static graphs, 
neglecting the temporal dimension [l6l - [l8j . The time 
evolution of a real system, when considered, is usually 
studied by evaluating the standard static measures (dis- 
tances and clustering coefficient) on snapshots of the net- 
work taken at different times [l9|, [2(| ■ As we will show 
below, this approach does not capture entirely the dy- 
namic correlations of a time-varying network. In this 
Letter, we introduce a measure of the temporal distance 
between the nodes of time-varying graphs, i.e. graphs 
with time-fluctuating links. Such temporal distance takes 
into account the actual time order, duration and corre- 
lations between links appearing at different times. This 
metric, together with a measure of the time-persistence 
of the links, allows to define and investigate temporal 
small-world behavior in social and biological networks 
that change over time. 



Consider a network with N nodes, where the links can 
fluctuate in time. The typical example is a social sys- 
tem, with no births or deaths, where the patterns of 
interaction are changing in time because of the spatial 
movement of the individuals, or because the individuals 
lose old acquaintances and get new ones. The system 
can be described at its maximum resolution sampling 
time as a time-varying graph, i.e. a discrete sequence 
(an ordered set) G\, G2, . . . , Gt of T undirected or di- 
rected graphs, where T is the length of the sequence. 
In compact notation, we denote the entire sequence as 
Q = {Gt}t=i,2,...,T- A time-varying graph Q can be rep- 
resented by means of a iV x N time-dependent adjacency 
matrix A(t), t — 1, . . . , T, where are the elements 

of the adjacency matrix of the i-th graph. We indicate 
as K(t) the number of links in the t-th graph of the se- 
quence. A sequence of graphs is convenient to describe 
systems where each connection starts at a specific time, 
and also has a temporal duration. In this sense, time- 
varying graphs are different from previous temporal ap- 
proaches [J [U, Hi], [22| designed to characterize systems 
as email exchanges, where the links have instead no tem- 
poral duration, because the exchange is instantaneous. 
Moreover, a time-varying graph is a different ensemble 
from those usually studied in the literature 0, [24j . In 
fact, in a time-varying graph, what matters is not only 
the probabilty distribution P(G) over the graphs in the 
ensemble, but also how the graphs arc ordered in time. 
By counting the number of times a given graph G ap- 
pears in the time sequence, we can construct P(G). To 
fully describe time-varying graphs we also need to know 
how graphs are correlated in time. For instance we need 
to know the conditional probabilities P(G t \G t -i) of ob- 
serving graph Gt after graph Gt-i (more in general, the 
probabilities P{G t \G\, G2, ■ ■ ■ , G t _i) of observing graph 
G t after the sequence G\, G%, . . . , Gt-i). In most cases, 
the contacts between the same node pair in time- varying 
systems tend to be clustered in time, i.e. they show per- 
sistence over time [J]. For instance, people tend to en- 
gage in relations for continuous intervals of time. Hence, 
a given link has a higher probability to appear in graph 
Gt if it was already present in graph Gt_i. To quantify 
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this effect, following Ref. [2(| we compute G, the average 
topological overlap of the neighbor set of a node between 
two successive graphs in the sequence: 

(i) 

We name this metric the temporal- correlation coefficient 
of (y. The value of C is in the range [0,1]. In particular, 
if all graphs in the sequence are equal, we have C = 1 . 

A fundamental concept in graph theory is that of 
geodesic, or shortest path. In a static graph, a shortest 
path between nodes i and j is defined as a path of min- 
imal length between the two nodes. This is a sequence 
of adjacent nodes starting at i, ending at j, and visiting 
the minimum number of nodes. Finally, the distance be- 
tween node i and node j is set equal to the length of the 
shortest paths from i to j. Here, we introduce the con- 
cepts of temporal shortest path and temporal distance to 
generalize the definitions of shortest paths and of node 
distance to the case of time-varying graphs. We illus- 
trate the basic idea with the example shown in Fig. QJi. 
Suppose node A wants to send a message in the fastest 
possible way to the other nodes of the graph. We assume 
that node A can start passing the message at time t = 1 , 
and the message has to be delivered by time t = 4. On 
graph Gi , node A can directly pass the message to nodes 
B and D, which are therefore assigned temporal distance 
1 from node A, since they can be reached in one unit of 
time. There are also other temporal paths to go from 
A to nodes B and D in three time units. For example, 
we can go from A to D in the following way: A — > B in 
G\, B — > D in G3. This is also a temporal path from 
A to D, though it is not the shortest, since the fastest 
way to go from A to D is to use the link A — > D in G\. 
Distance 3 is assigned to node C, since the message can 
be passed from A to D in graph G±, and then from node 
D to node C in G3, thus reaching G in three time units. 
Node F can be reached in 4 time steps by means of three 
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FIG. 1: An example of a time- varying graph Q with T = 4 
(panel a), and its projection into a static graph (panel b). 

alternative shortest paths: A — > B in G\, B — > F in G4; 
A-*-DwGi,D->C'wG 3 ,C->Fin G 4 ; and A -> D 
in G\, D — > B in G3, B — > F in G4. Finally, there are no 



temporal paths from A to E, hence we set the temporal 
distance of E from A equal to 00, and we say that E is 
not reachable from A. This is an effect of the time order 
of the links in a time-varying graph, and indeed node A 
and E are connected in the aggregate graph shown in 
Fig. [lb hr which all links arc considered as concurrent. 

Notice also that, due to the time order of the links, 
the temporal distances are not symmetric, even if the 
time-varying graph consists of a sequence of undirected 
graphs. For instance, while the temporal distance from 
A to F is 4, the temporal distance from F to A is 00 
(because the links occur in the wrong time order to facil- 
itate the passage from F to A). Conversely, in the static 
graph in Fig. [Tb, we have oIaf = dp a — 2. In order 
words, the main difference between a time- varying graph 
G, as that shown in panel a), and its associated static 
graph, reported in panel b), is that some of the short- 
est paths of the static graph are not temporally valid 
(in the sense that the links do not appear in the correct 
time order) and, therefore, cannot be used to route mes- 
sages. In general, in time- varying graphs there are more 
disconnected node pairs, than in static ones. As an ex- 
ample, the static graph in Fig. [T]is composed of a single 
connected component, while if time is taken into consid- 
eration, it is not possible to go from A to E, or from F to 
A. In order to compute the temporal distances dij for all 
node pairs i,j = 1, 2, N of a generic graph Q, we have 
implemented a generalization of the breadth first search 
algorithm. The average temporal connectivity properties 
of Q can be measured by the characteristic temporal path 
length L: 

Alternatively, in order to avoid the potential diver- 
gence due to pairs of nodes that are not temporally con- 
nected, we can define the temporal global efficiency of Q 
as E3: 

Low values of L (high values of E) indicate that the nodes 
of the graphs can communicate efficiently. In the follow- 
ing, we will show that time-varying graphs from models 
and real-world systems can be, at the same time, tem- 
porally clustered and still have small temporal distances 
between their nodes. In analogy with the small-world 
analysis in static graphs [HI, [13] , we will compare the ac- 
tual values of G, L and E of a given time- varying graph Q, 
with the corresponding values calculated by considering 
an ensemble {g rand } of randomized versions of Q. Each 
sequence g rand is obtained by randomly reshuffling the 
graphs in Q, i.e., by destroying the time order (and cor- 
relations) in the original sequence Gi, G2, . . . , Gt- More 
precisely, we will show that some time-varying graphs 
can have a value of G much larger than the correlation 
coefficient of the reshuffled sequence C rand , and, at the 
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same time a value of L as small as L rand . We will refer 
to this behavior as small-world behavior in time-varying 
systems. 



L(p.)/L(0) 



C(p.)/C(0) 
■ . j 



FIG. 2: Characteristic temporal path length and temporal- 
correlation coefficient of time- varying graphs produced by the 
model of moving agents, as a function of the probability pj of 
long-distance jumps. In the simulations we have set N = 100, 
D = 100 m, v — 1 m/s, r c = 5 m and produced sequences of 
length T = 500. The characteristic temporal path length of 
the reshuffled sequences is reported as dashed line. 

Random-walkers network model- We first illustrate 
how this behavior can be obtained in a network model 
of moving agents, as a result of simple motion rules. We 
consider a system of N random walkers which move in 
a two-dimensional square of linear size D with a fixed 
velocity v, and additionally perform long-distance jumps 
to randomly chosen position of the square with a jump 
probability Pj 3]. For each fixed value of Pj £ [0,1], the 
time-varying network Q is constructed by linking, every 
second, all nodes having a distance in space smaller than 
a given value r c . In Fig. [5] we plot C and L as a func- 
tion of pj. The values reported are normalized to the 
maximum values of C and L obtained for pj = 0, and 
respectively equal to C*(0) = 0.91 and L(0) = 442.8. We 
observe that a small percentage of jumps is sufficient to 
create links between nodes otherwise at large temporal 
distances and to produce a large drop in the temporal 
L. When pj = 0.01, L has reduced to one forth of L(0), 
and when pj — 0.1, L has about the same value as for 
the reshuffled sequence. The value of L rand obtained as 
an average over 1000 realizations of g rand is reported as 
dashed line. While L(pj) is rapidly decreasing, C(pj) is 
constant up to large values of pj ~ 0.1, so that for inter- 
mediate values of pj we have time- varying graphs exibit- 
ing small-world behavior. Finally, we have found that, 
by approximating for each value of pj the corresponding 
time- varying graph Q as a static graph, we obtain a value 
of static L (not reported in figure) which slowly changes 
in the interval [0, 1]. For instance, for pj = 1, we have 
L(pj)/L(0) = 0.74. We thus cannot capture the tempo- 
ral small-world behavior with the standard characteristic 
path length of a static graph. 

Brain cortical networks.- We finally explore real- world 
time- varying complex networks. We first consider time- 



varying functional cortical networks extracted from a set 
of high-resolution EEG recordings in a group of 5 normal 
subjects performing a task consisting in a foot movement 
[TTj ] . For each subject, and for each of four frequency 
bands (a, /?, 7, 9), we considered a time period of 0.5 sec 
corresponding to the final phase of execution of the foot 
movement. Each time-varying graph has N = 16 nodes, 
representing cortical regions of interest, and consists in 
a time sequence of T = 100 directed unweighted graphs, 
where the directed links represent causal influences be- 
tween cortical regions (see Ref. ll| for details). We have 
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0.000059 



TABLE I: Temporal-correlation, characteristic temporal path 
length and efficiency for brain cortical networks (subject 1, 
and four band frequencies) [Tl| . for the social interaction 
networks of INFOCOM'06 (time periods between 1pm and 
2:30pm, four different days), and for messages over Facebook 
online social network (three different months of year 2007) 
[2^ | . Results are compared with those obtained for 1000 ran- 
domized (shuffled) sequences of the same length. The values 
in parenthesis next to C ran d are the respective standard devi- 
ations. The values in parenthesis next to L and JJ and are the 
percentage of pairs of nodes that are temporally connected 
and not considered in the averages. 

computed the values of C, L and E for each real sequence 
and for the reshuffled ones. In Table Q] we report the re- 
sults for one of the subjects. For all the considered bands, 
the real sequence exhibits small- world properties, having 
a large value of C (significantly larger than C rand ) and, 
at the same time, a small characteristic temporal path 
length (a high efficiency), comparable to that observed 
in the shuffled sequence. Similar results (not reported) 
were obtained for the other four subjects. 

Social interaction networks.- The second real case 
study of our analysis is a time-varying social network 
based on a dataset of contacts among participants of 
INFOCOM'06, a major data communication conference 
which took place in a hotel. The contacts were collected 
by means of Bluetooth-enabled devices able to record in- 
teractions among people that are in promixity [26] . The 
discovery process of new devices was performed every 2 
minutes. In Table H] we report the data for the interac- 
tions during lunchtime between 1pm and 2:30pm. This 
is the interval with the larger number of contacts dur- 
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ing a day. Each sequence is made of T — 45 undirected 
unweighted graphs with N — 78 nodes each. The av- 
erage path length and the efficiency are similar for the 
original and reshuffled traces (the number in parenthesis 
close to L and L rand are the percentage of pair of nodes 
being temporally connected and hence considered in the 
computation of the average path length), whereas C is 
more than double that of C rand . This can be considered 
as an indication of small-world behavior in these traces 
according to our definition. 

Online social networks.- The third system we study is 
based on interactions over an online social network. The 
original dataset contains the messages sent among 6 mil- 
lions users in the London network of Facebook over one 
year (March 2007 to February 2008) We have di- 

vided the contacts according to the months of the year 
and, for each month, we have filtered out all contacts 
between pairs of nodes which exchange less than 10 mes- 
sages per month. This allows us to consider only the sub- 
set of most active users, obtaining networks with about 
N = 100, 000 users per month. For each month, the time 
varying graph is composed by T = 720 (or T = 744) di- 
rected graphs, one for each hour of the month. As shown 
in Table U for four different months of the dataset, the av- 
erage temporal path length of the time- varying networks 
is close to the value obtained for the reshuffled sequences. 
However, the network under study is disconnected in sev- 
eral different components, and only an extremely small 
percentage (about 10~ 6 ) of the node couples are tempo- 



rally connected. Consequently, the characteristic tem- 
poral path length was evaluated as an average over a 
small number of node couples. A better characterization 
of the system can be obtained by means of the tempo- 
ral efficiency. The values of E and E rand measured for 
Facebook are in general smaller than those observed in 
the other two networks, this being due to the high dis- 
connectedness of Facebook. Nevertheless, as for the case 
of the cortical networks and of INFOCOM'06, the real 
Facebook is almost as efficient as its reshuffled version. 
Finally, also for Facebook we observe a temporal small- 
world behavior: while the length of the temporal paths of 
the time-varying network are not affected by the reshuf- 
fling procedure, the temporal correlation coefficient C is 
about one order of magnitude larger than in the reshuf- 
fled version C rand . 

In conclusion, our results suggest that time- varying 
networks, strongly clustered in time and, at the same 
time, with short temporal paths between their nodes, 
might be widespread in biological, social and man-made 
systems, often with important dynamical consequences 
[25j . We hope that our work will stimulate further stud- 
ies of temporal small-world behavior in real time- varying 
systems. 
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